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ABSTRACT 

In this paper we present a numerical study of plasma jets produced by intense 
laser matter interactions. Through this study we hope to better understand 
astrophysical jets and their recent experimental simulations in the laboratory. 
We paid special attention to radiation cooling and the interaction of the jet with 
ambient gas. Four cases are presented in this paper; two of them deal with 
the propagation of jets in vacuum, while in the other two the propagation takes 
place in the ambient gas. Available experimental results are reproduced to good 
accuracy in the vacuum case. For jets in ambient gas, we find that the existence 
of the surrounding gas confines the jet into a narrow cylindrical shape so that 
both the density and temperature of the jet remain high enough for effective 
radiation cooling. As a result, a collimated plasma jet is formed in these cases. 
The dimensionless parameters characterizing the laboratory jets and protostellar 
jets have overlapping domains. We also discuss the cooling lengths for our model 
and compare them with the corresponding values in the astrophysical jets. A 
plasma jet in the ambient gas experiment is proposed which is within the reach 
of present day technology, and can be relevant to astrophysical phenomena. 

Subject headings: ISM: jets and outflows — hydrodynamics — plasmas — meth- 
ods:numerical 



1. Introduction 

Astrophysical jets found in supersonic bipolar flows from black holes, active galactic 
nuclei(AGN), protostars and similar phenomena are subject to active research. Several 
problems concerning jets are remaining to be solved, for example, why do the jets propagate 
such a long distance so stably. Many numerical simulations of astrophysical jets have been 
performed so far. 2D non-relativistic calculations have been done for AGN jets, by Norman 
et al. (1982). They found some dimensionless parameters relevant to the propagation of jets 
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over a long distance. 2D and 3D relativistic calculations have been carried out by Marti et 
al. (1997); Aloy et al. (1999). Protostellar jets got more attention, detailed numerical simula- 
tions have been performed by several authors (see Reipurth, Bally & Devine (1997); Reipurth 
et al. (1997); Reipurth & Hearthcote (1997); Bally & Reipurth (2001) for observational data). 
Since the radiation cooling and/or external magnetic field are supposed to have very impor- 
tant effect on the plasma jets propagation, radiation hydrodynamic/magnetohydrodynamic 
codes are essential. For example, Blondin et al. (1990); Stone & Norman (1993) showed 
with their radiation hydrodynamic codes that the radiation cooling plays an important role 
for the collimation, leading to the formation of a dense cold shell at the head of the jet. 
Todo et al. (1993) studied, by means of with MHD code the propagation of jets from young 
stellar objects into ambient matter, and showed that the azimuthal magnetic field of ~ 70 
fiG is strong enough to drive a helical kink instability. Since then, many efforts have been 
put to simulations in 2D and 3D with both radiation cooling and magnetic fields included 
(Cerqueira & Gouveia Dal Pino 1997; Frank et al. 1998; Gardiner et al. 2000; Frank et al. 
2000). The parametric researches that vary the form of the injection, steady or pulsed jets, 
and the configuration of magnetic field (helical or longitudal) helped to clarify the systemat- 
ics of jets. In particular, it was shown that a helical magnetic field can affect the morphology 
of the head of jets. It was also demonstrated that Kelvin-Helmholtz instability is crucial for 
the morphology of jets, and linear as well as nonlinear stability analyses were carried out 
with or without magnetic field to characterize the influence of these factors on the jet shape 
(Xu et al. 2000; Hardee k Stone 1997; Stone & Hardee 2000; Ryu et al. 2000). 

So far these studies are purely theoretical. Recently, however, with the advance of fem- 
tosecond laser technology, an experimental investigation of these astrophysical phenomena in 
the laboratory seems to be possible (Remington et al. 1999, 2000b). This new scientific field 
is called laboratory astrophysics or laser astrophysics. This progress enable the investigation 
of radiation hydrodynamics, hydrodynamic instabilities, and similar astrophysical phenom- 
ena in model experiments with intense lasers. Such experiments might be able to give us a 
new and improved insight into those phenomena which could not be obtained by studies of 
extraterrestrial observations or numerical simulations. Some experiments along these lines 
have already been carried out, in which hydrodynamic instabilities and the interactions of 
a strong shock with matter have been measured. Other experiments have been proposed 
(Remington et al. 2000a). Backup by numerical simulations are indispensable, because they 
can connect between the processes occurring in laser produced plasmas to astrophysical phe- 
nomena that are the aim of our studies. This paper is one of those attempts along this line. 
In this paper, we will propose a different experiment to understand astrophysical jets. 

Some experiments relevant to astrophysical jets have been performed with intense lasers 
(Farley et al. 1999; Shigemori et al. 2000; Stone & Turner 2000). The so-called 'cone' target 
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(Fig.l) is employed in those experiments. Another type of experiment has been recently 
done by Lebedev et al. (2001), using a conical array of fine metric wires. In general, intense 
laser matter interaction generates high temperature and high density plasma which starts to 
ablate immediately. If the intensity of laser is ~ 10 14 W cm -2 , the expansion velocity of the 
plasma reaches a few hundreds of kilometers per second. This is of the same order as the 
velocity of protostellar jets. This indicates that jets which realistically mimic astrophysical 
ones can be produced in laboratory. In these experiments, however, the plasma expands into 
vacuum, and one can not the investigate influence of ambient gas, which does exist for the 
astrophysical jets. As a matter of fact, the ratio of the mass density of the jet to that of the 
interstellar gas (p j e t I P ambient) is 1 ~ 10 for protostellar jets. 

In this paper, we investigate the propagation of laser-produced jets not only in the 
vacuum but also in the ambient gas using the 2D hydrodynamics code developed recently 
(Mizuta 2001). Although the code has two versions, relativistic and non-relativistic ones, in 
the present paper we use only the non-relativistic code. 

This paper is organized as follows. In sections 2 and 3, the basic formulas, the numerical 
method and the initial conditions are presented. The main results are shown in section 4. 
We discuss some implications to astrophysical jets in section 5 and give a short summary in 
section 6. 



2. Basic Equations and Numerical Method 

Assuming a plasma having axial symmetry 2D Euler equation for perfect fluid is solved: 

du ld(rf) dg ^ 

dt r dr dz 

where the conserved variable vector u, the flux vectors / and g, and the source vector s are 
defined as, 



u= (p,pv r ,pv z ,E) T , (2) 

/ = (Pv r , pvl + P, pV r V z , (E + p)v r ) T ', (3) 

9 = (pv z , pv r v z , pv 2 z +p,(E + p)v z ) T , (4) 

s = (0,p/r,0,0) T , (5) 

E = p(e + (v 2 + v 2 z )/2). (6) 



In the above equations, p, Vi,p, and e are mass density, i-component of velocity, pressure, and 
specific internal energy, respectively. We ignore viscosity and heat conduction of the fluid. 
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These assumptions are justified as follows. As for the collisions, the parameter 5 = \ m f p / L sys 
defined as the ratio of the collisional mean free path X m f p of electrons or ions and, the 
characteristic length of the system L sys , is 10~ 4 ~ 10~ 6 for our models described below. 
Hence a hydrodynamic description of the jets is appropriate. Other dimensionless parameters 
characterizing the dissipative processes are the Reynolds number Re = TjVj/v which is a 
measure of the local viscosity and the Peclet number Pe = r-jVj/uh which describes the heat 
conduction inside the jet. Here rj, Vj, v and v h are the radius and velocity of the jet, the 
viscosity and thermal diffusivity coefficients respectively. Criteria for the influence of these 
dissipative phenomena on the flow are Re < 1 and Pe < 1. In our models we find Re > 10 9 
and Pe > 10 5 . Thus the neglection of these dissipative processes is well justified. These 
numbers are evaluated at t — 2.0 ns at the center of the inflow region, see below. 

An ideal gas equation is used in the computations, 

P = P<7 - 1), (7) 
where 7 is the adiabatic exponent (=5/3 in this paper). 

We adopt the Marquina's flux formula (Donat & Marquina 1996) which is based on an 
approximate Riemann solver derived from the spectral decomposition of the Jacobian matrix 
of Euler equations. We have developed a computational code which solves Eq. (1) and (7) 
The code uses cylindrical coordinates with 300(r) x 1500(z) grid points. The accuracy of 
the code is the second order in space due to a MUSCL method(van Leer 1977, 1979) and is 
the first order in time. 

It is assumed that the plasma is optically thin, so the radiation cooling effect is impor- 
tant. In the cooling term | J ra d\, only for the bremsstrahlung is included (Zel'dovich & Raizer 
1966), 

I J rad \ = 1.42 • IQ- 27 Z* 2 T l ' 2 n e n ion eig cnrV 1 , (8) 

where Z*,T,n e and n ion are the average ionic charge, temperature and number densities of 
electron and ion, respectively. The temperature and average ionic charge are related to the 
specific internal energy e by, 

3_(z* + i)t 

C "2 m t ' [) 

where m; is the mass of ion. The average ionic charge state Z* and temperature T were cal- 
culated iteratively from local instantaneous density and specific internal energy, as obtained 
from the hydrodynamic calculation by using Eq.(9) and approximate expression based on the 
Thomas Fermi model for Z* (Salzmann 1998). The formulas for the line emission processes 
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were computed in the usual manner and are not included here for simplicity. As already 
mentioned the radiation transport is also neglected in this paper. Since we are interested 
mainly in the qualitative effects of radiation cooling, these approximations are adequate. 

We have to distinguish the target matter from the ambient matter in order to make 
the radiation cooling effective only for the former. For this purpose we introduced another 
continuity equation for the target matter: 

W) | 1 d(rpv r f) | djpvj) _ Q 
dt r dr dz 
where / is the fraction function. This equation is solved simultaneously with the hydro- 
dynamic equations. With Z*,T, and / thus obtained, the cooling term is calculated in a 
separate step. For numerical convenience, the radiation cooling term is turned off for very 
low temperatures, less than 70 eV. The radiation cooling term is neglected also in the vicinity 
of the target, where the radiative processes are very complicated and the basic picture of 
our model is not valid anyway. 



3. Initial Condition 

In order to simulate the experiments, we put a 'cone' target at an end of the compu- 
tational region. At t — 0, the laser deposits all its energy on the surface. The depth of the 
surface is ~ 10/im. A Gaussian shaped laser pulse of duration of fOOps was assumed. This is 
much shorter than the dynamical time scale, which is more than a few nanoseconds. When 
the target made of gold is irradiated by intense laser of ~ 10 14 W cm -2 , the temperature of 
the target rises to a few hundred eV up to I keV, and the ablation takes place. Because the 
velocity of the ablated plasma is at most its sound velocity, v a b ~ c s ~ 10 7 cm s _1 , the ab- 
lated plasma expands only up to ~ 10/zm during the laser irradiation ~ f OOps. This is much 
smaller than the scale of the system we consider in this paper. The laser energy deposited 
into thermal energy is E\ = 526J for all simulations. Energy of ionization is not included, 
therefore in a real experiment the incident laser energy, of course, has to be significantly 
larger. 

Four simulations have been done in this paper. In cases 1 and 2, the ambient gas is very 
dilute (p a = I0~ 6 g cm -3 ). In the first case, we neglect the radiation cooling, while in the 
second case the radiation cooling is turned on. These intend to simulate the experiments. In 
the other two cases, 3 and 4, we let the plasma expansion in the presence of an ambient gas. 
Case 3 neglects the radiation cooling, and case 4 is done with the radiation cooling term. 
No corresponding experiments have been carried out as yet. These models serve to clarify 
the effect of the ambient gas. 
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The temperature of the target surface, where the laser energy is deposited, increase to 
~ 400 eV. The angle of the cone is 126°, which is similar to the experiments (see Fig. 2). 
The target matter is chosen to be gold ( 79 Au) so that we can obtain large effect of the 
radiation cooling due to a high-Z plasma. The density of the target p t is initially uniform, 
p t = 19. 2g cm -3 , and equals to the solid density of the gold. Table 1 summarizes the 
parameters of the initial conditions. 



4. Results 
4.1. Vacuum case 

First, we discuss cases 1 and 2 which simulate plasma expansion into vacuum. This 
condition is similar to the experiments of Farley et al. (1999); Shigemori et al. (2000). For 
numerical reasons, we assumed in these cases as well a cold (T a = 0.026eV) and very low 
density (p a = 10~ 6 g cm~ 3 ) ambient gas. Due to this very low density this ambient gas had 
no influence on any of our results. Figure 3 shows the density contours at different times 
with the radiation cooling off (Fig. 3a) and on (Fig. 3b), respectively. 

In the case of no radiation cooling (Fig. 3a), we can see two regions, that is, the inflow 
region in which matter flows to the axis, and the outflow region. This is schematically 
shown in Fig.4. Most of the laser-heated matter expands from the target, forming plasma 
with velocity of a few hundreds of kilometers per second, farther away from the target it 
converges to the symmetry axis, making a nozzle-like structure in the inflow region. Then, 
the plasma in the nozzle spouts out from the tip of the nozzle, forming the outflow region. 
Because the ambient matter is very dilute, there is no pressure support to sustain this nozzle 
structure. As a result, the outflow spreads out in all directions. This flow structure is suitable 
for low Z targets, because the effect of the radiation cooling is indeed negligible. 

An analysis of the nozzle structure in the inflow region is in order. The converging 
flow from the cone target to the central axis is very similar with the structure considered 
by Canto et al. (1988). From their theoretical model we can estimate the beam radius r,- 
and velocity Vj from the converging flow angle 9, the reciprocal of the compression ratio £ 
between the density of the nozzle to that of the inflow and the velocity of the converging flow 
v (see Fig. 5). These are 9 = 26.5° and £ = 1/4 in our case. Then r,- and Vj are determined 
as 

tana 

r i = I — V* U ) 

tan 9 + tan a 

cos(# + a) 

vj = v , (12) 

cos a 
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where yo is the width of the inflow and a is the angle between the conical shock and z-axis 
and given as 

(l-Q + Kl-^-^tan 2 ^/ 2 

tana = . 13 

2tan# v ; 

Applying this theory to our case with y = 250/iin and v o ~ 3 x 10 7 km s~\ we obtain 
Tj = 180/im and Vj ~ 1 x 10 7 km s _1 , in good agreement with our numerical results. 

For high Z targets, we have to consider radiation cooling effects. When radiation cooling 
is taken into account, the flow structure is changed dramatically (Fig.3b). The plasma in 
the inflow region is collimated strongly and shows a very thin nozzle at a later time (the 
radius of the structure is ~ 40/xm). Although the collimation of the plasma is sustained in 
the inflow region by converging inflow from the target, the matter again spreads out in the 
outflow region, resulting in almost the same structure as that for the case without radiation 
cooling. These experiments are not suitable for the study of the propagation of collimated 
jets. 

Figures 6 a, b and c are the density, pressure and temperature profiles along r-axis, 
respectively, at z — 1000/im, t = 2.0ns. In the case with radiation cooling, the plasma is 
cooled efficiently around the symmetry axis, leading to the increase of the pressure gradient. 
As the density increases, the plasma is further cooled, and as a result a jet-like nozzle 
is formed, which was actually observed in the experiments. Altogether, our simulations 
support the experimental results that the radiation cooling is a crucial ingredient for the 
formation of this nozzle. 



4.2. Dense gas case 

In cases 3 and 4, the ambient matter density was much higher, p a = 10~ 3 g cm -3 . The 
time evolution of hydrodynamics leading to the formation of the inflow and outflow regions 
is essentially unchanged. Figure 7 shows the density contours at different times with the 
radiation cooling off (Fig. 7a) and on (Fig. 7b), respectively. The bow shock, which accelerates 
the ambient matter, can be seen clearly. 

In the case with radiation cooling, the flow structure is dramatically changed as is 
obvious in Fig. 7b. The effect of radiation cooling manifests itself in both the inflow and the 
outflow regions. In the inflow region, the nozzle structure becomes very thin just like in the 
vacuum case in the later time. The main difference appears in the outflow region. Unlike 
all the other cases, the width of the outflow region is very small. Opposed to intuition, this 
narrow structure is not a direct outcome of the thin nozzle in the inflow region. This is 
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understood from the second panel of Fig. 7b, in which there is already a very thin outflow 
region while there is no collimated structure in the inflow region. When the matter enters 
the outflow region, it tends to spread out as in vacuum case. However, the shocked ambient 
matter slows the expansion in the direction perpendicular to the symmetry axis, and the 
density and temperature of the flow remain high enough for the radiation cooling to be 
effective. As a result, a collimated plasma is produced in the outflow region. 

This difference of structure in the outer region is remarkable, as is further evident in 
Fig. 8, in which the density, pressure and Mach number profiles along r-axis, respectively, at 
z = 2000/im, t = 4.0ns are shown. The radius of jet is about 40/zm, and the density ratio is 
V = Pjet/Pa ~ 100, namely, this is a so-called dense jet. From Fig. 8b, we see that the narrow 
structure in the outflow region is sustained by the shocked ambient matter, whose pressure 
is about 1000 times larger than that of the unshocked ambient matter. The jet is supersonic 
Mj et ~ 70. The velocity of jet is ~ 400km s -1 . This is of the same order as the velocity of 
protostellar jets. In order to study the propagation of astrophysical jets in the laboratory, 
we think it is necessary to do experiments with ambient gas. 

5. Discussion 

Firstly we consider the scaling law between a laboratory jet and a protostellar jet. The 
velocity of these two jets is almost the same-a few hundred kilometers per second. The ratio 
of the time scales is, however, widely different I0 2 ~ 4 yr for protostellar jets and 10 _9 s for 
laboratory jets, a scale up of a factor of 10 18-20 . For the experiment to simulate efficiently 
the astrophysical jet, the ratio of the lengths of jets has to be ~ I0 18 ~ 20 . The typical ratio 
of the radii of the jets (Vp^/r^) is I0 15 cm/10~ 3 cm = I0 18 . Thus, the aspect ratios of the 
jets are close to each other. The temperature is a few eV for both jets, which is consistent 
with the scaling of velocity mentioned above. Thus, a good scaling law holds between the 
protstellar jet and the laboratory jet considered in this paper. 

The relevant parameters which characterize the jet are the density ratio rj, the Mach 
number and velocity of jets Mj et , Vj et , the pressure ratio K and the cooling parameter x- As 
the magnetic fields were not included our code, no comparison can be made with magneto- 
hydrodynamic simulations (Todo et al. 1993; Cerqueira & Gouveia Dal Pino 1997; Frank et 
al. 1998; Gardiner et al. 2000; Frank et al. 2000). As is obvious from the results obtained 
both for the vacuum and the ambient gas cases, the radiation cooling is the most crucial 
ingredient in the collimation of the jets. 

The cooling parameter is defined from the cooling term which has already been employed 
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in the hydrodynanmic simulations, \J ra d\ in Eq.(8). It is the ratio 

cooling length v jet T rad 

X = — T-7 T- ~ , 14 

jet radius Vj et 

where Vj et and rj et are the velocity and radius of jet, respectively, and r rad is a characteristic 
radiation cooling time. r rad is defined from the ratio of thermal energy density to emission 
power as 

Trad = T^—y 
I -Jrad] 

The smaller x is> the more effective the cooling is. The cooling term adopted in Eq.(8) 
is different from those used in other papers (Blondin et al. 1990; Stone & Norman 1993). 
We assumed in this paper that the jet matter is Au, rather than hydrogen in most other 
papers. Our choice was influenced by the need to enhance the radiation cooling so that the 
hydrodynamics is affected within a few nanoseconds. In fact, the average ionic charge state 
Z* is about 20 — 40 for T = 50 — lOOeV. Although the detailed dependence on T is different 
between our cooling term and others, it is emphasized that the dynamics should be similar 
as long as the cooling parameter x has a similar value during the greater part of the time 
evolution. 

In case 4 (with dense ambient gas and radiation cooling), the typical radiation cooling 
time is T rad ~ 10ns at t — 1.0ns, and r rad ~ 100ns for t = 4.0ns. These were evaluated at the 
points with the largest emission power. Accordingly, the radiation cooling parameters for 
each time are x ~ 10 and x ~ 100, respectively. Thus, the radiation cooling is effective for 
the earlier phase, while the jet is cool enough, while the radiation cooling is ineffective for 
the later phase. However, it is noted that we underestimate the radiation cooling, since we 
neglected the line emissions which could be substantial in our cases, and take into account 
only the bremsstrahlung. This was done mainly to keep numerical simplicity. 

The aim of our calculations and future experiments is to understand better the physics 
of protostellar jets. Although the observational values of physical quantities are rather 
uncertain, they are typically, n e ~ 10 3 cm -3 , Tj ~ 10 4 K, Vj ~ 100km s _1 and r 3 - = 10 15 cm 
(Bally & Reipurth 2001). These values roughly correspond to rj ~ 1 — 10 3 ,M_ ?ei ps 1 with 
the main uncertainty coming from the lack of data on the density and temperature of the 
ambient matter. From these parameters, it is estimated that the total emission power is 
\Jrad\ ~ 10~ 18 erg cm~ 3 s -1 for collisional by excited hydrogen atoms (Dalgarno & McRay 
1972). The emission power due to the bremsstrahlung is | J rad \ ~ 10~ 19 erg cm _3 s _1 . Thus 
the bremsstrahlung is substantially weaker, and in general the bound-bound emission is 
dominant in such plasmas. The emission power | J rad \ ~ 10~ 18 erg cm~ 3 s _1 gives a cooling 
time of r rad ~ 10 9 s, which then leads to a cooling parameter of x ~ 10. It is evident from 
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these values that the experimental values used in this paper are appropriate for studying 
astrophysical jets. As mentioned in the introduction, Blondin et al. (1990); Stone & Norman 
(1993) studied the properties of propagating jets by varying r\ between 1 and 10. Their 
parameters correspond to x between 0.1 and 10. We think that our results have overlapping 
region to the initial conditions of those papers. This means that the experiments proposed 
here might give some new insights into the astrophysical jets studied theoretically in their 
papers. 

6. Conclusion 

In this paper results of four hydrodynamic simulations have been shown, in which jet 
were generated from 'cone' target irradiated by intense laser. In all the simulations, we have 
found a common feature, that is, the existence of an inflow region and an outflow region. 
Depending on whether the radiation cooling is turned on or off, and whether there is some 
ambient gas or not, the flow structure inside these regions is very different from each other. 

When the plasma flows into a very low density ambient gas, the collimation of the plasma 
jet occurs only due to radiation cooling. With radiation cooling, the collimated plasma shows 
up in the inflow region. However, since the structure is sustained by converging inflow from 
the target, it disappears in the outflow region, in which there is no pressure support to 
maintain it. This feature is very similar to that found in the experiments. We have also 
found that the nozzle structure with no radiation cooling is well described with the analytical 
model by Canto et al. (1988) 

In the case with ambient gas, a bow shock appears very clearly. The main difference 
between the vacuum and dense gas cases is that in the latter case the collimated flow appears 
not only in the inflow region but also in the outflow region. Particularly, the collimated 
plasma in the outflow region is sustained by the pressure of the shocked ambient matter. 

We have shown the possibility to carry out a laboratory experiment, using high intensity 
laser-plasma interaction, to generate collimated plasma jets propagating into an ambient 
gas. The parameters of such plasma jet are similar to those of protostellar jets. If such an 
experiment is performed, we can get important information about the physics of astrophysical 
jets. 

External magnetic fields were not included in our calculations. Experiments with mag- 
netic fields have already proposed and will be carried out in the future. We intend to adjust 
our simulations to the needs of future experiments. 
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Table 1: Detailed parameters for initial condition, p and T mean density and temperature. 
Index V and 'a' stand for target and ambient matter 





easel case2 case3 case4 


target 


An 


ptg cm~ 3 


19.2 


T t eV(hot/cold) 


676/1.77 


p a g cm" 3 


1(T 6 io- 3 


T a [eV] 


0.026 


cooling term 


OFF ON OFF ON 
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Fig. 1. — 'Cone' target. (a):Some intense lasers irradiate on the 'cone' target surface, 
(b): Ablation plasma will collimate and become a jet-like structure. 
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Fig. 2. — Initial conditions. The cone target is set at the left end of computational region 
(1200/xm x 3000/iin). Hot plasma is put at the surface. It's depth is ~ 10/im. 
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Fig. 3. — Density profiles in r — z plane, in the vacuum case; (a)radiation 
off (easel ),(b)radiation on(case2). t =0.0(1), 1.0(2),2.0(3),2.6(4),3.0(5)ns 
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Fig. 4. — Schematic flow structure. The flow is divided to two regions, inflow and outflow 
region. 
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5. — A schematic drawing of the conical flow to the central axis. 
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Fig. 6. — Density, pressure and temperature profiles along r axis at z=1000/xm, t=2.0ns. 
The radiation cooling term is OFF(dash line,case 1) and ON(solid line, case 2). 
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Fig. 7. — Density profiles in r — z plane, in the dense g ; (a)radiation 

off(case3),(b)radiation on(case4). t =0.0(1), 1.0(2),2.0(3),3.0(4),4.0(5)ns 
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Fig. 8. — Density, pressure and Mach number profiles along r axis at z=2000 /im, t=4.0 ns. 
The radiation cooling term is OFF(dash line, case 3) and ON(solid line, case 4). 



